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With the identification of (_D3.j(2317), I>sj(2460)) as the (0+, 1+) doublet in the heavy quark 
effective field theory, we derive the light cone QCD sum rule for the coupling of eta meson with 
Dsj{2'in)Ds and Dsj(2460)-Ds . Through r; — ■k'^ mixing we calculate their pionic decay widths, 
which are consistent with the experimental values (or upper limits). Combining the radiative decay 
widths derived by Colangelo, Fazio and Ozpineci in the same framework, we conclude that the decay 
patterns of Dsj (2317, 2460) strongly support their interpretation as ordinary cs mesons. 
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I. INTRODUCTION 



, Since the discovery of D=7(2317) H and D.^.y ( 2460) 0, there have been lots of experimental investigations of these 

iiiiiiiiiimiii 

e observed final states. T 

prediction which are really unexpected. Many theoretical papers have been dedicated to the understanding of 



two narrow resonances |MIJilMIS.SMHIiyili3- They have the natural spin-parity assignment as the 0+, I'*' charm- 
strange mesons from the observed final states. Their masses are about one hundred MeV lower than the quark model 



>' their underlying structure [11 Hill III [11 [11 [ii 113, 111 E 

0^ Proposed schemes include the (O"'",!"'") chiral partners of the D.,, Dt doublet in heavy quark effective field theory 
^ ' [13|, lilJllLpfs: molecules [3^ fourquark states [ll[ll[ll[2l|2ll[2l[2l[2|,i:)7r atom ^ and conventional cs 

Q ! states[lMl[S[2l[2l[Mlla[sl[lllH3 

' These two states are lower than the DK and D*K thresholds respectively. Their strong decays are isospin violating 
[ and occur through two steps: Dsj{2317) ^ Ds + v ^ D, + Tr° and Dsj(2460) Dl +_n-^ D* + 7r°. The second 
■ step is induced by the rj — tt^ mixing due to the mass difference between to,, and 1371. There have been some 
discussions of their strong and radiative decays within the quark model [isL IT^ Is^ Is^ 133. l4d. l4]| . The strong decay 
widths from various approaches differ significantly as can be seen from Table ^ Their radiative decay widths were 
r^' calculated using light cone QCD sum rule (LCQSR) not long ago 42]. Very recently, the branch ratios of the st rong 
O-i: and radiative decays were measured quite accurately by Belle Collaboration [1 Q and Babar Collaboration [1 [lTI| ! 
In order to pin down the underlying quark content of these narrow states, a reliable calculation of their strong decay 
widths will be very helpful. 

In this work, we assume (Dsj(2317), iDsj(2460)) as the cs states and study their strong decays in the LCQSR 
framework, which has been used extensively in extracting low-lying hadron masses and coupling constants in the past 
decade (see Ref. ^3] for a review). This paper is organized as follows. We calculate the coupling constant QD^aDs-q 
d ' and the strong decay width of I?sj(2317) Dgir^ through r/ — 7r° mixing in Section^ _Dsj(2460) decay is presented 
in Section IIIII We compare our results with experimental data and other theoretical approaches in literature and 
summarize our results in Section HVI . We collect the light cone wave functions of the r] meson in the appendix. 



II. Dsj(22,n) ^ Ds ^ r) ^ Ds + -k" 
The amplitude of the strong decay L's,/(2317) Dg + rj can be defined as 

iV (q) Ds{p)\D so {p + q)) mD,o9D,oD,r, ■ (1) 

We calculate the coupling constant gDsoD^ri through the following correlation function 

F{p^{p + qf)^^ I d^x eP-%n{q)\T[4{x)UQ)m . (2) 
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where Jo (a;) = c{x)s{x) and J^ix) = c{x)i^5s{x) are the interpolating currents of Dgo and Dg respectively. 

At the quark level, the correlation function can be expressed in terms of the eta meson light cone wave functions 
after the insertion of the charm quark propagator at the leading order 



(0|T{c(x)c(0)}|0)=*^°(a;) = 



(2n)H 



I? — fc^ — ie 



Now we have 



(2^ 



d^x '—^r^ Ue(r;(g)|s(x)i75s(0)|0) - ifc"(?7(g)|s(a;)7„s(0)|0) 



m 



(3) 



(4) 



In order to include the contribution from the twist-four eta meson light-cone wave functions, we need the three 
particle piece in the charm quark propagator: 



(0|T{c(x)c(0)}|0) = 



d^k 
(2^ 



dv 



1 (i + mc 
2(m2-A:2)2 



G^''ivx)a^, + 



1 



ll - fc2 



vx^G^^vx)-i, 



,(5) 



where G^i, = G'^^^ with tr(A°A'') = 26'^^, and gg is the strong coupling constant. The complete expression of 
F(p2^ + g)2-| twist-four reads: 



du 



ml-{p + uq)' 



r| - Fjjiprf{u)p.q + rricFri^ijjippiu) 



-F,,m2(l H — £ - 

' ''^ ml-{p + uq) 



:)G(u) + -F„ml — - — — tttII , „ , ,„ , 



A{u)p.q 



2mt 



+Fr,ml I dv I Va 



r7i2 — (p + ?i(7)2 m2 — (p + uqy 

2_(^^(^.^ ,.^.\n\2 {M'^vi- - n) - (2'/'-L - n + - ^ii)} . (6) 



p.q 



' nil ^ iP+ (o^i + '"Ct3)qy 



where 



G{u) = duB{u) . 



(7) 



friiu), ipp{u), B(u) etc are light cone amplitudes of 77 meson |4J,|45|, which arc collected in the Appendix, i^,, — ~-^f, 
with frj defined as 
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(u(0)7;.75?/(0) + d(0)7^75d(0) - 2s(0)7^75s(0)) ^(q)) = if^q^ 



At the phenomenological level, F{p^, [p + qY) can be expressed as 



F{p\{p + qf) 



(8) 



(9) 



The ellipse denotes the contribution from the continuum. The decay constants Jdso ^'^d fo, are defined as 

{0\jI\Ds) = Id mo, , (10) 



(11) 



Applying the double Borel transformation with respect to p^ and [p + q^ to Eqs. l|njl and © and invoking the 
quark-hadron duality, we get the following sum rule: 



foMfo^gD^nD^T, 



-e 2i\ 



,2 



iM [e ^-e ^] -Af F^(p^{uo) + mcFni-i,,ipp{uo) 



F^m''^G{uo) - ^F^mlA'iuo) + ^F,,ml{2h(2^j_ - - Wv± - n + 2<^± - "^n)) 



-e 'm^ Fjfm'^m'^ 



G{uo) + {uo) 



(12) 



«o = l/2 
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The functions Ii and I2 are defined as 

1 



Uq 



dai\ J^{ai,l-uo,uo~ ai)- da^ ^ 

Uo-ai Juo-ai "3 

i-uo jr(uQ,l - uo - a3,a3) , r° , F{ai,l ~ uq,uo - ai) 

da^ h / dai , 

"3 Jo "0 - "1 



(13) 
(14) 



Since there exists an overlapping working window for the two Borel 



where J- is one of the twist-four hght cone amplitudes Lp\\, ip±, i^y, ip±. 

In Eq.Ca, A'P = ^0. , uo = . 
parameters Mf,M|, it's convenient to let = M|, i.e., uq — i. The eta meson light-cone wave functions are 
known quite well at mq = 5. Such a choice allows the clean subtraction of the continuum contribution. We can simply 



introduce a threshold parameter so and replace M^e with 



to subtract the contribution from 



the continuum and excited states |4g . 

In the numerical analysis, the values of parameters in the above sum rule are: mn ^r, — 2.317 GeV, m/)^ = 1.968 
GeV, nic (1 GeV) = 1.35 GeV, fo, = (266 ± 32) MeV 113, /d,o = (225 ± 25) MeV For fo, and fo^o wc use the 
central values. Values of the other parameters are given in the Appendix. The variation of ^Dso-Dsj; with M'^ for the 
different sq is shown in Fig. ^ 
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FIG. 1: The variation of the coupling constant go^oDsv with (in unit of GeV'^). The long-dashed, short-dashed and solid 
curves correspond to so — 6.0, 6.25, 6.5 GeV'^ respectively. 



In the working region of the Borel parameter 5 GeV^ < < 7 GeV^, we get 

-3.02 < go^aD^,^ < -2.66 , 



(15) 



where the uncertainty arises from the variation of Af^ and sq. Numerically, the twist-three term ipp has the largest 
contribution to the sum rule. 

The pionic decay of Dsj(2317) occurs through 77 — tt*^ mixing, which is described by the isospin violating piece in 
the chiral lagrangian 



2 j-2 



r 



"4/ 



4(m„ -f mrf) 



Tr(em,e + ^V,e^) 



(16) 



where ^ = exp(j7f//^), tt the light meson octet and is the light quark mass matrix. The mass difference between 
up and down quarks induces the 77 — tt" mixing with a suppression factor around — . Finally the strong decay 

ii^— i 

width reads 



r(i?,,;(2317) ^ i?.^°) = ^.g2,^^^^^( 
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Numerically we have 



T{Ds J (2317) ~* DsTT^) = (34 - 44) keV 



(17) 



(18) 
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III. J (2460) ^Dt+v^ D*s + 
For Dsj(2AQ0) decay, we define the following matrix element 

{l]{q)D*{p)\D',i{p + q)) = ™i5^^5d:,D*^'?^ 



(19) 



where rj^^ and are the polarization tensors for the 1+ and 1 states D^i^D* respectively. We start from the 
correlation function 

F,.{p\{p + qf) = z I d^xe^'^-''{r,{q)\T[jlix)J^m\0) (20) 
where Jfi{x) — c{x)jfjs{x) and J^{x) — c{x)^ulbs{x)- At the hadron level, we have 



FMP^ip + qV) 



vyp- ' 2m^, ml. 



■qt^Pu + • • •) 



(m|)^ -p2)("^|,^„ - {p + qYY 

where we have kept the g^u and q^Pu structures. The decay constants and Jbz are defined as 

{0\J+\D:) = fD-mn^^e^ , 
(Z?li|J,^|0) = f^,m^,y,. 

Following the same procedure as in Section^ we obtain a sum rule from the g^i, structure 



^2 1^2 



-e 2M^ J M fe — e "m^ 



+ ^Fr,mlA{uo) + i^^^m2(2/i(^||) - hi^^^ + 2^^ - ^||))] 



+e "iw^ 



tio = l/2 



Similarly we can get a second sum rule from the q^Pu structure 



2 1^2 



^D'jDt9D:,Dtr, = —2 Zt^^ ^ i [e ^ - 6 ^] (mq)] 



+e "iij^ 



1 1 777- 

-mcF,^firjip„{uo) + -to^F^(1 + -^)A(mo) - 2F^rn^MoG(uo) 



+F^m2(jg(^ll ^ 2(^_L - ^||) - 2/4(2(/j_l + '^||)) 



110 = 1/2 



The functions /s and I4 in Eqs. (|24|l and (|25(l are defined as 



da I 



da I 



^ T{ai,l - ai - a^.as) 
das , 

^ ~ "1 1 ^ 

das 2 — -^("1; 1 - ai - ct3, as) . 



(21) 



(22) 
(23) 



(24) 



(25) 

(26) 
(27) 



Unfortunately the sum rule Eq. (|25|l is very unstable. There is no working window for the Borel parameter M^. 
In the following we focus on the sum rule Eq. We use m^'^ = 2.460 GeV, to_dj = 2.112 GeV jig, fjj'^ ^ 

fD^ — fo, E3- The variation of gjy' ^jy.^j with is presented in Fig. 2. 

In the working window of Borel parameter 4 GeV^ < < 6 GeV^ , we have 

2.61 < 5^^^^.,, < 3.14 . (28) 
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FIG. 2: The variation of the coupling constant 



with (in unit of GeV'^). The long-dashed, short-dashed and soUd 



curves correspond to so — 6.25, 6.5, 6.75 GeV respectively. 



The contribution from ipp term is also very important numerically. The pinonic decay width reads 



Finally we have 



1447r 



-(2 



)(- 



2 



r(L>^j(2460) ^ + tt") = (35 - 51) keV 



(29) 



(30) 



IV. DISCUSSION 

The strong decay widths of _Ds,7(2317) — > Dg + tt'^ and Dsj{2460) D* + 7r° have been calculated by several 
groups. Their results are collected in Table ^ together with ours. The first five calculations assume cs picture while 
the last two use composite non-cs pictures. The decay width of _Ds,/(2317) is roughly the same as that of Z3s,/(2460) 
from all approaches. The 1/toc correction is expected to modify the small values in the second column from vector 
dominance model in the heavy quark limit j43 | . 

TABLE I: Strong decay widths (in keV) of Dsj (2317) and 73s j (2460) from various theoretical approaches. 







LCQSR 




[13] 


[38j 


[39] 


[17] 


m 


0:.7(2317) - 


->73s7r" 


34-44 


7± 1 


21.5 


~ 10 


16 


10-100 


150 ± 70 


J (2460) - 


- 7?>° 


35-51 


7± 1 


21.5 


~ 10 


32 




150 ± 70 



The radiative decay widths of 73^,7(2317,2460) were calculated using LCQSR in [4a|: r(735j(2317) D* + -f) = 
(4 - 6) keV, r(L>,j(2460) ^ 73^7) = (19 - 29) keV, r(73,j(2460) ^ 73^ + 7) = (0.6 - 1.1) keV, r(73^j(2460) ^ 
Dsj{2317) + 7) = (0.5 - 0.8) keV. Experimentally only 73sj(2460) 73^7 has been observed by Belle H |3| and 
Babar H, 0|. We have collected the experimental ratio of radiative and strong decays of Dsj mesons together with 
the central values of theoretical predictions from LCQSR based on Ref. 42] and present work in Table ^ For 
73^7(2460) Dgj), we get a range 0.37-0.83 for the ratio, consistent with both Belle and Babar's measurement. 

In short summary, we have calculated the coupling constants gDsoD^ri and Qj^i ^.^^ in the framework of LCQSR. 

Through the ry - tt" mixing we obtain r(735,7(2317) 73,7r°) ^ (34 - 44) keV^and r(735,7(2460) DItt°) = 
(35 — 51) keV. These two widths are similar in magnitude, as expected from heavy quark symmetry. The ratio between 
the radiative widths and strong decay widths obtained in the same LCQSR framework is consistent with Belle and 
Babar's most recent data, which strongly indicates 73^7(2317) and 73^7(2460) are conventional cs mesons. In the future, 
B decays into 73^7 mesons may also play an important role in exploring these charming states [igllsollsllfs^lsMls^lssll . 
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TABLE II: Comparison between experimental ratio of Dsj (2317, 2460) radiative and strong decay widths and theoretical 

I ' 

predictions from LCQSR based on Ref. ^42j and this work. 





Belle 


Babar 


CLEO [2] 


LCQSR 


r(D*,(2317)^D;7) 
r(l3J^(2317)^Os7r") 


< 0.18 [4] 




< 0.059 


0.13 


r(i3sj(2460)^I3s7) 
r(D3j(2460)->DJirO) 


0.55 ±0.13 ±0.08 [4] 


0.375 ± 0.054 ± 0.057 [11] 


< 0.49 


0.56 


r(l33j(2460)^D'7) 
r(D3j(2460)^D;7rO) 


< 0.31 \£ 




< 0.16 


0.02 


r(Ds J (2460)^D*j (2317)7) 




< 0.23 |1D1 


< 0.58 


0.015 


r(D3j(2460)^Dj7rO) 
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Appendix 



We use qTq to denote {uTu + dTd— 2sTs) /V^- Up to twist four, the two- and three-particle light-cone wave functions 
of eta meson can be written as 113, Ha : 



< ??k(2;)7/. 759(0) |o > 



^ifvQf^ J du e™«^ [if,^ {u) ± —mi^x^A{u)] 



2 ■ 'qx Jo 
< v\Qix)i-f5qiO)\0 > = fn^^n / du e''"^^ifp{u) , 



< ri\q{x)af,^j5q{0)\0 > = -frifJ-niqf.x^ - q^x^,) J due ™'>''ip„{u), 

< v\q{x)<^al3759sG^,^{ux)q{x)\0 > = ^f,,^J.rJV3[iq^iqa9ly|3 - qvqa9tj,ti) - {qt^qpgiya ~ quqi39tj.a)] 



<v\q{^)l^l^.gsG^p{vx)qmO> = /^m2[g^(.g„^-^) j Va^v^{a,)e~''^<''^+-"^'^ 



q ■ X / . 



+Uml^{q^xp-q0X^) / Pa,(^|| (aOe"*'"'"^+""'^ 



< r]\q{x)j^gsGap{vx)q{x)\0 > = j/,,m^ 



2 qi±_ 

X 



qp\9o 



'^aq^, 

q ■ X 



qa I gpf. 



q ■ X 



ifrjm^n — iqcxp - q0Xa) / I?a,^||(aOe^^"("^+'^"^^ , 
'q-x J 

where the operator Gaf3 is the dual of Gap'- Gap — ^^a/sSpG^''; Vai is defined as Vai = daida2daz5{l — ai — a^ — as) 
and ~ 1.2/^ = 0.156 GeV, % = 0.013, to^(1 GeV) = 0.125 GeV, m,, = 0.548 GeV, ^ 2.4 GeV 113. 
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The distribution amplitudes ipri etc can be parameterized as 

(t>p{u 

&u{l - u){\ + (5% - ^773^3 - - lpla2)]CY\0 , 



9r,{u 



6uu{l + aaCf '(C) + 04^^ '(C)) , 
5 
2 



1 + (30773 - y,)C',^'{C) + him^3 - - ^pf,a2)Cy'(C) 

1 27 n 3 



</'3r,(ai, ^2, as) 
V9||(ai, 02,013) 

((Cj^(ai, 02, as) 
(^11 (01,02, as) 
<f± (ai,a2,03) 



^ _A6 24 20 

6WW{— + + 2O773 + y?74 

18 

+ ( a2 + 21r/4a)4){2M^(10 - 15m + 6m') Inu + 2i2^(10 - 15u + 6u') Inu 

5 

+Mi2(2 + liuu)} , 
360?73aiO2O3{l + ^^3(703 - 3)} , 
120020103(010(01 — 02) , 

3003(02 - oi)[/ioo + hoiaz + ^/iio(5a3 - 3) , 
120010203(^00 + Wio(3o3 - 1) , 

3 

-SOogl/iooOs + /ioi[oga3 - 601O2] + ft.io[o3a3 - -(o^ + al)] , 



where u=l — M,C = 27i— l,o = l — o. C2 4''^^' 

gives the mass correction and are defined as = Oy, Vij and hij are related to hadronic matrix elements 774, W4 
and 02 as 



(C) are Gegenbauer polynomials. Here gri{u) = B{u) + (pri{u). 
ij , Vij and hij 



aw = ^774^4 - Jja2 , vw = ^m^i 



, Woo = -3774, 
^01 = ^774^^4 - 4a2 , ft-io = i?74'^4 + 4a2, voa = -hvi- 



The values of 02 et al are: 02 = 0.115, 04 = —0.015, 773 = 0.013, uj^ = —3, 774 = 0.5, 0^4 = 0.2. All of them are scaled 
at ;u = 1 GeV. 



[1] Babar Collaboration, B. Aubert et al., Phys. Rev. Lett. 90, 242001 (2003). 

[2] CLEO Collaboration, D. Besson et al., Phys. Rev. D, 032002 (2003). 

[3] Belle Collaboration, P. Krokovny et al., Phys. Rev. Lett. 91, 262002 (2003). 

[4] Belle Collaboration, Y.Mikami et al., Phys. Rev. Lett. 92, 012002 (2004). 

[5] Belle Collaboration, A. Drutskoy et al., Phys. Rev. Lett. 96, 061802 (2005). 

[6] Belle Collaboration, P. Krokovny et al., hep-ex/0310053 

[7] FOCUS CoUabortation, E. W. Vaandering, hep-ex/0406044 

[8] Babar Collaboration, B. Aubert et al., Phys. Rev. Lett. 93, 181801 (2004). 

[9] Babar Collaboration, B. Aubert et al., Phys. Rev. D. 69, 031101 (2004). 
[10] Babar Collaboration, G. Calderini et al., hep-ex/0405081 
[11] Babar Collaboration, B. Aubert et al., hep-ex/0408067 

[12] S. Godfrey and N. Isgur, Phys. Rev. D 32, 189 (1985); S. Godfrey and R. Kokoski, Phys. Rev. D 43, 1679 (1991); M. Di 

Pierro and E. Eichten, Phys. Rev. D 64, 114004 (2001). 
[13] W. A. Bardeen, E. J. Eichten and C. T. HiU, Phys. Rev. D 68, 054024 (2003). 
[14] M. A. Nowak, M. Rho and L Zahed, Acta Phys. Polon. B 35, 2377 (2004). 
[15] E. Kolomeitsev and M. Lutz, Phys. Lett. B 582, 39 (2004). 
[16] T. Barnes, F. E. Close, H. J. Lipkin, Phys. Rev. D 68, 054006 (2003). 
[17] H.-Y. Cheng and W.-S. Hon, Phys. Lett. B 566, 193 (2003). 
[18] K. Terasaki, Phys. Rev. D 68, 011501 (2003). 

[19] T. E. Browder, S. Pakvasa, and A. A. Petrov, Phys. Lett. B 578, 365 (2004). 



8 



U. Dmitrasinovic, Phys. Rev. D 70, 096011 (2004); Phys. Rev. Lett. 94, 162002 (2005). 

A. Hayashigaki and K. Terasaki, h ep-ph/0411285, 

M.E. Bracco et al., hep-ph/050313'^j 

H. Kim and Y. Oh, hep-ph/050825l' 

M. Nielsen et al., hep-ph/0509131 

A.P. Szczepaniak, Phys. Lett. B 567, 23 (2003). 

A. Deandrea, G. Nardulh, A. D. Polosa, Phys. Rev. D 68, 097501 (2003). 

R. N. Cahn and J. D. Jackson, Phys. Rev. D 68, 037502 (2003). 

Y. B. Dai, C. S. Huang, C. Liu and S. L. Zhu, Phys. Rev. D 68, 114011 (2003). 

W. Lucha and F. Schobert, Mod. Phys. Lett. A 18, 2837 (2003). 

J. Hofmann and M.F.M. Lutz, Nucl. Phys. A 733, 142 (2004). 

M. Sadzikowski, Phys. Lett. B 579, 39 (2004). 

D. Becirevic et al., Phys. Lett B 599, 59 (2004). 

L W. Lee et al., hep-ph/0412210 

A. Zhang, Phys. Rev. D 72, 017902 (2005). 

UKQCD Colaboration, A. DongaU et al, Phys. Lett. B 569, 41 (2003). 

P. Colangelo, F. De Fazio and R. Ferrandes, Mod. Phys. Lett. A 19, 2083 (2004). 

P. L. Cho and M. B. Wise, Phys. Rev. D 49, 6228(1994). 

S. Godfrey, Phys. Lett. B 568, 254 (2003). 

Fayyazuddin and Riazuddin, Phys. Rev. D 69, 114008 (2004). 

T. Mehen and R. P. Springer, Phys. Rev. D 70, 074014 (2004). 

P. Colangelo and F. De. Fazio, Phys. Lett. B 570, 180(2003). 

P. Colangelo, F. De Fazio and A. Ozpineci, hep-ph/0505195 

P. Colangelo and A. Khodjamirian, in 'At the Frontier of Particle Physics/Handbook of QCD', ed. by M. Shifman (World 

Scientific, Singapore, 2001), hep-ph/0010175j 

P. Ball, JHEP 9901, 010 (1999) 

P. Ball, R. Zwicky, Phys. Rev. D 69, 115001 (2005) 

V. M. Belyaev, V. M. Braun, A. Khodjamirian, R. Riickl, Phys. Rev. D 51, 6177 (1995). 
S. Eidelman et al. [Particle Dada Group], Phys. Lett. B 592, 1 (2004). 
S. Ishida et al, hep-ph/0310061 
M. Suzuki, hep-ph/0307118 

C. H. Chen and H. N. Li, Phys. Rev. D 69, 054002 (2004). 
H. Y. Cheng, Phys. Rev. D 68, 094005 (2003). 
A. Datta and P. J. ODonnell, Phys. Lett. B 572, 164 (2003). 
M. Q. Huang, Phys. Rev. D 69, 114015 (2004). 

H. Y. Cheng, C. K. Chua and C. W. Hwang, Phys. Rev. D 69, 074025 (2004). 
T. Barnes, hep-ph/0406327| 



